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Abstract 



Spherically symmetric expansionfree distributions are systematically stud- 
ied. The whole set of field equations and junction conditions are presented 
for a general distribution of dissipative anisotropic fluid (principal stresses 
unequal), and the expansionfree condition is integrated. In order to under- 
stand the physical meaning of expansionfree motion, two different definitions 
for the radial velocity of a fluid element are discussed. It is shown that 
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the appearance of a cavity is inevitable in the expansionfree evolution. The 
nondissipative case is considered in detail and the Skripkin model is recov- 
ered. 
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1 Introduction 



The problem of general relativistic gravitational collapse of massive stars has 
attracted the attention of researchers for many years, starting with the semi- 
nal paper by Oppenheimer and Snyder [1]. The motivation for such interest is 
easily understood: the gravitational collapse of massive stars represents one 
of the few observable phenomena, where general relativity is expected to play 
a relevant role. Ever since that work, much was written by researchers try- 
ing to provide models of evolving gravitating spheres (see [2] and references 
therein). However this endeavour proved to be difficult. Different kinds of 
obstacles appearing, depending on the approach adopted for the modelling. 

Thus, numerical methods allow for considering more realistic equations 
of state, but the obtained results, in general, are restricted and highly model 
dependent. Also, specific difficulties, associated to numerical solutions of 
partial differential equations in presence of shocks, complicate further the 
problem. Therefore it seems useful to consider nonstatic models which are 
relatively simple to analyze but still contain some of the essential features of 
a realistic situation. For doing so we need to appeal to a simple equation of 
state and/or to additional physically meaningful heuristic assumptions. In 
this work, we shall assume the fluid to be expansionfree. 

As is well known, the motion of a fluid may be characterized by the four 
acceleration vector (a a ), the shear tensor (<7 aj g), the expansion scalar (0) and 
the vorticity tensor (which vanishes in the spherically symmetric case). The 
relevance of the shear tensor in the evolution of selfgravitating systems and 
the consequences emerging from its vanishing has been brought out by many 
authors (see [3] and references therein). 

In this work we shall study the properties of an expansionfree spherically 
symmetric selfgravitating fluid. 

Since the expansion scalar describes the rate of change of small volumes 
of the fluid, it is intuitively clear that the evolution of an expansionfree 
spherically symmetric distribution should necessarily imply the formation of 
a vacuum cavity within the distribution (see a more rigorous argument on 
this in section V). Thus, in the case of an overall expansion, the increase in 
volume due to the increasing area of the external boundary surface must be 
compensated with the increase of the area of the internal boundary surface 
(delimiting the cavity) in order to keep G vanishing. The argument in the 
case of collapse is similar. 
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For sake of generality we shall start our discussion by considering an 
anisotropic dissipative viscous fluid (arguments to justify such kind of fluid 
distributions may found in [4]- [6] and references therein). For this kind of 
distribution we shall write the field equations, the junction conditions, at, 
both, the inner and the external boundary surface (section II), as well as 
the dynamical equations (section III). Next we shall integrate the expansion- 
free condition and find the general form of the metric for the anisotropic 
dissipative viscous fluid (section IV). 

In order to understand better the physical meaning of the expansionfree 
motion, we shall discuss two different definitions of radial velocity of a fluid 
element, in terms of which both the expansion and the shear can be expressed 
(section V). 

We shall next consider the nondissipative case, and we shall specialize 
further to the isotropic fluid. In this latter case we shall recover as a particular 
example the Skripkin model [7], assuming the energy density to be constant 
(section VI). 

Finally a discussion on the results is presented in the last section. 

2 The energy-momentum tensor, the field equa- 
tions and the junction conditions 

We consider a spherically symmetric distribution of collapsing fluid, bounded 
by a spherical surface Y^ e \ The fluid is assumed to be locally anisotropic 
(principal stresses unequal) and undergoing dissipation in the form of heat 
flow (to model dissipation in the diffusion approximation), null radiation (to 
model dissipation in the free streaming approximation) and shearing viscos- 
ity. 

Choosing comoving coordinates inside T^ e \ the general interior metric 
can be written 

dsl = -A 2 dt 2 + B 2 dr 2 + R 2 (d6 2 + sin 2 6d<f) 2 ), (1) 

where A, B and R are functions of t and r and are assumed positive. We 
number the coordinates x° = t, x 1 = r, x 2 = 9 and x 3 = <fi. Observe that A 
and B are dimensionless, whereas R has the same dimension as r. 
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The matter energy- momentum T a/3 inside X* 6 ) has the form 

T~p = (v+P±)V a V /3 +P ± g a p+(P r -P ± )xaXp+<laVf3+V a qp+6lJp-2r)a a p, (2) 

where // is the energy density, P r the radial pressure, P± the tangential 
pressure, q a the heat flux, e the energy density of the null fluid describing 
dissipation in the free streaming approximation, rj the coefficient of shear 
viscosity, V a the four velocity of the fluid, x a a un it four vector along the 
radial direction and l a a radial null four vector. These quantities satisfy 

V a V a = -l, V a g a = 0, x a x a = h x a v a = o, l a V a = -l, l a l a = 0. (3) 

The acceleration a a and the expansion 6 of the fluid are given by 

a a = V a . if) Vf > , G = V a ;a . (4) 

and its shear a af3 by 

<r a p = V"( Q;/3 ) + a^Vp) - ^9{g a p + V a Vp), (5) 

We do not explicitly add bulk viscosity to the system because it can be 
absorbed into the radial and tangential pressures, P r and P±, of the collapsing 
fluid [8]. 

Since we assumed the metric (1) comoving then 

V a = A~ l 5^ q a = qB' l 5^ l a = A~ 1 5q + B- l 5*, X a = B^d^ (6) 

where q is a function of t and r. 

From (4) with (6) we have for the acceleration and its scalar a, 

A' ( A' \ 2 

°i = ^4' a2 = a,a(la= [ab) ' ^ 

and for the expansion 

where the prime stands for r differentiation and the dot stands for differen- 
tiation with respect to t. With (6) we obtain for the shear (5) its non zero 
components 

(Tii = ^B a, a 22 = -r-%- = --R a, (9) 
3 sm 9 3 



and its scalar 
where 



v H cr a(3 = -a , 



}_(B_R 
A\B~R, 



2.1 The Einstein equations 

Einstein's field equations for the interior spacetime (1) are given 



and its non zero components with (1), (2) and (6) become 
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The component (14) can be rewritten with (8) and (10) as 
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Next, the mass function m(t,r) introduced by Misner and Sharp [9] (see 
also [10]) reads 



R 3 

- ft d 23 

m = —R23 = — 



R 
~2 




(18) 



Thus in the most general case (locally anisotropic and dissipative) we have 
available four field equations (13-16) for eight variables, namely A, B, R, //, 
P r , P ± , e and q. Since we are going to consider expansionfree systems we have 
the additional condition O = 0. Evidently, in order to find specific models 
(to close the system of equations) we need to provide additional information, 
which could be given in the form of constitutive equations for q and e, and 
equations of state for both pressures. 



2.2 The exterior spacetime and junction conditions 

Outside £ (e) we assume we have the Vaidya spacetime (i.e. we assume all 
outgoing radiation is massless), described by 



ds 2 = - 



1 - 



2M(v) 



dv 2 - 2drdv + r 2 (dd 2 + sin 2 6d(f) 2 ), 



(19) 



where M(v) denotes the total mass, and v is the retarded time. 

The matching of the full nonadiabatic sphere (including viscosity) to the 
Vaidya spacetime, on the surface r = r E ( e ) = constant, was discussed in [11]. 
From the continuity of the first and second differential forms it follows (see 
[11] for details) 



m(t,r) S = M(v), 



(20) 



and 



IK BR' _RA" 
\~R ~ ~ ~R~A , 
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] R 



R\ R 
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' A_ R!\ R^ 

2 H + ~r) ~R 



R) 



21 



(21) 



where = means that both sides of the equation are evaluated on T,^' (ob- 
serve a misprint in eq.(40) in [11] and a slight difference in notation). 
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Comparing (21) with (14) and (15) one obtains 

We) 4 

q = P r - - V a, (22) 

Thus the matching of (1) and (19) on E (e) implies (20) and (22), which 
reduces to equation (41) in [11] with the appropriate change in notation. 
Observe a misprint in equation (27) in [5] (the a appearing there is the one 
defined in [11], which is —1/3 of the one used here and in [5]). 

As we mentioned in the introduction, the expansionfree models present 
an internal vacuum cavity. If we call E^ the boundary surface between the 
cavity and the fluid, then the matching of the Minkowski spacetime within 
the cavity to the fluid distribution, implies 

m(£,r) S =0, (23) 

vM 4 

q = P r - 3 w- (24) 



3 Dynamical equations 

To study the dynamical properties of the system, let us introduce, following 
Misner and Sharp [9], the proper time derivative Dt given by 

^41- < 25 > 

and the proper radial derivative Dr, 

D *=hl- < 26) 

where R defines the areal radius of a spherical surface inside Y,^ ( as mea- 
sured from its area). 

Using (25) we can define the velocity U of the collapsing fluid (for another 
definition of velocity see section 6) as the variation of the areal radius with 
respect to proper time, i.e. 

U = D T R < (in the case of collapse). (27) 
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Then (18) can be rewritten as 

^ B! 
E = — = 
B 

With (26) we can express (17) as 
4vr(g + e) = E 



1 + U 2 - 



2m y/ 2 
R J 



Using (13)-(15) with (25) and (26) we obtain from (18) 

4 



Dtitl 



-4tt 



(p r + e- -rja^j U + (q + e)E 



R 2 , 
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DrUI = 4tt 
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which implies 
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47T 
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R dR 



(28) 



(29) 



(30) 



(31) 



(32) 



(assuming a regular centre to the distribution, so m(0) = 0). 

Expression (30) describes the rate of variation of the total energy inside a 
surface of areal radius R. On the right hand side of (30), (P r + e — 4r]cr/3)U 
(in the case of collapse U < 0) increases the energy inside R through the 
rate of work being done by the "effective" radial pressure P r — 4r]cr/3 and 
the radiation pressure e. Clearly here the heat flux q does not appear since 
there is no pressure associated with the diffusion process. The second term 
—(q + e)E is the matter energy leaving the spherical surface. 

Equation (31) shows how the total energy enclosed varies between neigh- 
boring spherical surfaces inside the fluid distribution. The first term on the 
right hand side of (31), /i + e, is due to the energy density of the fluid ele- 
ment plus the energy density of the null fluid describing dissipation in the 
free streaming approximation. The second term, (q + e)U / E is negative (in 
the case of collapse) and measures the outflow of heat and radiation. 

The non trivial components of the Bianchi identities, T7^ = 0, from (12) 
yield 



rp-a/3-rr _ 



{H + e) +(ji + P r + 2e- 



-r/a 
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+2^i + P ± + e + 2 -r)G 
(AR)> 



R 
R 



(g + e)' + 2(g + e)- 



(q + e) + 2(q + e) ( ^ + ^ 



AP 

B ' R 



= 0, 



(33) 



+ 



B 



(P r + e- ^qaj + (v + P r + 2e- ^qa 



+2(P r -P ± + e-2r)a)^- 



~A 

= 0, 



(34) 



or, by using (7), (8), (25), (26) and (28), they become, respectively, 



D T (fi + e) + -(3/2 + P r + 2P ± + 4e) 6 
o 



+-(P r -P ± + € - 2r)o)o + ED R (q + e) + 2(q + e) (a + 



E 
R 



0, (35) 



D T (q + e) + -(g + e)(20 + a) + ££> fl (p r + e - -77a) 
+ (/1 + P r + 2e - |t7ct^ a + 2(P r - P ± + e - 2^)^ = 0. (36) 

This last equation may be further tranformed as follows: The acceleration 
D T U of an infalling particle inside £ can be obtained by using (7), (15), (18) 
and (28), producing 



m 



D T U = -— -Au[P r + e- - W j R + Ea, 
and then, substituting a from (37) into (36), we obtain 



(37) 



/i + P r + 2e - -r]a ) D T U 



= -(fi + P r + 2e-- W 



-E 1 



m 



^—+4TT[P r + e--r]a)R 



D R (P r + e--r)a)+2(P r -P ± + e-2r)a)- 



-E 



Djiq + 2(7/ (2^ • 



(38) 
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The physical meaning of different terms in (38) has been discussed in detail 
in [4]-[6]. Suffice to say in this point that the first term on the right hand 
side describes the gravitational force term. 



4 Shearing expansionfree motion 

If the fluid has no expansion, i.e. 6 = 0, then from (8) we have 



B _ R 
B ~ ~ 2 i?' 



or, by integrating 



B = 



9(r) 
R 2 ' 



(39) 



(40) 



where g(r) is an arbitrary function of r. 
Substituting (39) into (14) we obtain 



R' n R R' R A' , , . „ 
R +2 RR-RA =i ^ + ^ AB ' 



(41) 



which can be integrated for R ^ producing 

/R 
(q + e)AB^dr 
R, 



R 2 R 
A = exp 



(42) 



where Ti(t) is an arbitrary function of t. With (40) and 42) then (1) becomes 

2 



ds 2 = - 



R 2 R 



cxp 



/R 
(q + e)AB^dr 
R 



dt 2 



+ 



— V 

R 2 ) 



dr 2 + R 2 (d9 2 + sin 2 6d(f) 2 ), 



(43) 



which is the general metric for a shearing expansionfree anistropic dissipative 
fluid. 

In the nondissipative case q = e = we can write (43) as 

dg2 = _ /^R\ df2 + l_ dr2 + R2 ^ 02 + §in2 0d(f)2 ^ (44) 
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where without loss of generality (by reparametrizing r) we put g — 1 (observe 
that a unit constant with dimensions [r 4 ] is assumed to multiply dr 2 ). Then 
we have that (44) is the general metric for a spherically symmetric anisotropic 
perfect fluid undergoing shearing and expansionfree evolution (observe that 
it has the same form as for the isotropic fluid [12]). 



5 On the physical meaning of expansionfree 
motion 

We shall now analyze under which conditions a time dependent spherically 
symmetric configuration may evolve without expansion. For doing so we shall 
try to develop more our understanding of the different speeds that involve the 
description of expansion as well as shear in the evolution of a selfgravitating 
fluid. The following discussion heavily relies on the kinematic quantities 
characterizing the motion of a medium presented in [13], with slight changes 
in notation. 

In Gaussian coordinates, the position of each particle may be given as 

x a = x a (y a ,s), (45) 

where s is the proper time along the world line of the particle, and y a (with 
a running from 1 to 3) is the position of the particle on a three-dimensional 
hypersurface (say E). Then for the unit vector tangent to the world line (the 
four-velocity) we have 

V = (46) 

and observe that 

! = (47) 

Next, for an infinitesimal variation of the world line we have 

5x* = —5y*, (48) 

from which it follows 

D T (5x a ) = V. a p 5x p . (49) 
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Introducing the projector hp on £ by 

hf = S$ + V a V 0: (50) 

we can define the position vector of the particle y a + Sy a relative to the 
particle ?/" on E, as 

5 ±x a = h%5x p . (51) 
Then the relative velocity between these two particles, is 

u a = h^D T (5 ± x p ), (52) 

and considering (49) and (51) it follows that 

u a = V. a p 5 L x p . (53) 

Now, the infinitesimal distance between two neighboring points on S is 

SI 2 = g a pS±x P S±x a , (54) 

then 

SlD T (Sl) = g a p5 ± xPD T (5 ± x a ), (55) 
or, by using (49) and (52), 

5lD T (5l) = V a . p 5 ± x p 5 ± x a . (56) 

Then, taking into consideration the expression for the irreducible components 
of a timelike vector 

V a -/3 = a a p - a a Vp + ^Qh a(3 , (57) 
where we assumed zero rotation, and substituting into (56) we obtain 

5lD T (5l) = 5 ± x^5 ± x a [a aP + , (58) 

or, introducing the spacelike unit vector 

13 



it becomes 



Dt(SI) _ „ , 



S] e a e"a a(3 + -. (60) 

Let us now consider, the spherically symmetric case, and apply (60) to 
two neighbouring points along the radial direction. In this case we have 
e a = x a , and using (6), (9) and (10) in (60) we obtain 



or, by using (11) and (8) 



51 3 3' 

D T (5l) B 



51 AB' 
Then with (27) and (62) we can write 



(62) 



D T (5l) D T R D T (5l) U 

a = —i b- = ^t-r> (63) 

and 

_ D T (6l) 2D T R _ D T (6l) 2U 

51 R 51 R' { ' 

Thus we see that in general there are two different contributions to the shear 
(63) and to the expansion (64). One is due to the "circumferential" velocity 
U [14], which is related to the change of areal radius R of a layer of matter, 
whereas the other is related to Dt(51), which has also the meaning of "ve- 
locity" , being the relative velocity between neighboring layers of matter, and 
can be in general different from U. 

From (63) we see that, if the spherical distribution of matter is collapsing, 
U < 0, the shear can vanish only if the relative distance between different 
layers of matter diminishes, Dt(51) < 0, and cancels the circumferential 
velocity. 

From (64) we see that the evolution of the fluid will be expansionfree, 
whenever the "circumferential" term cancels the term related to the vari- 
ation of distance of neighboring particles. Thus the collapse will proceed 
expansionfree, if the decrease of the perimeter of a comoving sphere (U < 0) 
is compensated by an increase in the distance of neighbouring particles (along 
the radial direction) according to (64). Alternatively, if the fluid is moving 
outward (U > 0) neghbouring particles will get closer (D T (5l) < 0). These 



14 



observations clarify further the origin of the cavity in expansionfree models. 
Indeed, consider two concentric fluid shells in the neighborhood of the cen- 
tre. As it follows from (74), close to the centre we have U ~ R. Now, in 
the process of expansion (increasing of R), the = condition implies as 
mentioned before that Dt{61) < 0, i.e. both shells become closer, however 
this would not be so as long as U ~ R, implying thereby that the O = 
condition requires that the innermost shell of fluid should be away from the 
centre, initiating therefrom the formation of the cavity. 

Let us see this from another perspective. Consider the infinitesimal vol- 
ume of the shell between two concentric spheres of radii r and r + 5r, 

5V = AnBR 2 5r, (65) 

then it follows 

D T (5V) = 47i(D T B)R 2 5r + 87rBR(D T R)5r, (66) 

or, dividing (66) by 5V and using (62), 

D T (6V) D T (6l) U 
5V ~ 51 + R- 

Which of course coincide with (64) since we know that D T (5V)/5V is the 
definition of the expansion. 

Thus in the process of contraction (expansion), the elementary volume 
5V decreases (increases) by two factors. On the one hand by the decreasing 
(increasing) of the areal radius R and on the other hand by the decreasing 
(increasing) of the proper radial distance between the two concentric sur- 
faces. Again, we have an expansionfree evolution (the elementary volume 5V 
remains constant) whenever the two contributions on the right hand of (67) 
cancel each other, in spite of the fact that neither of them vanishes. 

We shall now see how these two different definitions of radial velocity 
considered above are related. 

Let us first assume that U — 0. Then from (63) and (64) it follows O = a, 
feeding this back into (29) we get at once 

D T (U) = -^§±1>«, (68) 
thus, U — implies D T (5l) = only in the dissipationless case, q = e = 0. 
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Next, let us assume Dt{61) = 0. Then it follows from (63) and (64) that 
= —2(7, feeding this back into (29) we get 



D R a 



a 
R 



4vr(g + e) 



E 



whose integration with respect to R yields 

C 47T 



(69) 



(70) 



where ( is independent of R (for any layer of fluid, characterized by r = 
constant, ( may depend on t. In general however it may depend on t and r). 
On the other hand, (70) with (63) implies 



U = -( + 4tt J (q + e 



R 

E^ 



dR. 



(71) 



Since U — > as R — > 0, we must put ( = (if only the centre of the fluid 
distribution is covered by the coordinate system). Thus, in the nondissipative 
case U — 0. 

Therefore only in the nondissipative case U — <-> D T (5l) = (with the 
condition mentioned above). 
Finally, we can write (29) as 



U 



R 



D R [-) = —(q + e) + - 



E 



R 



which after integration with respect to R becomes 

cR 



or, 



U 



U = £R + R f 
Jo 



47T a 

¥<« + f) + S 



dR, 



R 



-R-R 



(e) 



47r . cy 
-E (q + C) + R 



dR. 



(72) 



(73) 



(74) 



In the shearfree nondissipative case we have from (74) that U ^ R, which 
is characteristic of the homologous evolution [15]. This implies that for two 
concentric shells of areal radii Ri and R 2 , we have in this case 



Ri 

Ro 



= constant. 



(75) 
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The second term on the right of (74) describes how the shear and dissipation 
deviate the evolution from the homologous regime. 

It is worth noticing that in the shearfree nondissipative case the sign of 
U for any fluid element is the same as that of £/ S ( e ) . However if we relax any 
of those conditions, that might not be true. Thus it would be possible, for 
example, to have an expandig outer shell with an imploding inner core. Such 
possibility was brought out before, but restricted to the quasistatic regime 
[16]. Here we see that such an scenario is also possible in the general dynamic 
regime. 



6 Shearing expansionfree perfect fluid 

We shall now restrict our study to a shearing expansionfree fluid without 
dissipation, q = e = and r\ = 0, then the metric reduces to (44) and the 
field equations (13-16), using (29), become 

1 T 2 

87r/i = -2R 3 R" - 5R 2 R' 2 + — - 3-^, (76) 

R l Ry 

l DRa+ R~ = > (77) 

8vrP r = -X ( 3* - 2*) + R 3 R> ( 2^ + 5^1 - -1 (78) 
R 5 R \ R n \ R R I R 2 



R 5 R V R n) 

while the Bianchi identities (33) and (34) read, 



R" R R R" , n R'\ 

+ 7^— + 3— + 10 — 
R R R R \R 



(79) 



/i + 2(PL-P r )| = 0> (80) 

P; + (fi + P r )S- + 2(/i + 2P r -P±)^ = 0. (81) 
R rl 

From (80) we have that if the fluid is isotropic, P r = Pj_, then the energy 
density fi is only r dependent. 

We can now integrate (76) under the assumption \i = [J,(r), to obtain 

2 _ 1 r 2 -2m rj 2 
R + + W (82) 
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where T2(t) is an arbitrary function of t and (32) has been used. 

We shall now specialize further our model to the case of constant energy 
density [7] 



6.1 The Skripkin model 

In [7] it is not explicitly assumed that 6 = 0, instead it is assumed that the 
fluid is nondissipative, has its energy density \i = /i = constant and the 
pressure isotropic. Of course, these conditions imply, because of (35), that 
= 0. Thus we have only one physical variable (P r ) and the system of field 
equations is closed and can be integrated. 

From the condition (j, — /i — constant, (82) becomes 

i?,2 = ~^ + ^ + § + ^' (83) 

with 

(84) 

It should be observed that (83) imposes a maximum to the value of R 
(Rmax), for which R' — 0. The physical origin of this maximum for the areal 
radius may be explained as follows: 

In the Skripkin picture the fluid is initially at rest, then there is a sudden 
explosion at the centre producing the outward ejection of the fluid, always 
keeping the conditions of nondissipation, ji = fj, = constant and isotropic 
pressure (i.e. = 0). Under these conditions, (38) becomes 

(/i + P) D T U = - (/i + P) + AtiPr}) - E 2 D R P, (85) 

with P r = P± = P ■ Now, as the fluid moves outward and R approaches 
Rmax, R' tends to zero which implies, because of (28), that E approaches 
zero. Furthermore, this implies that the gravitational (negative) term in (85) 
will prevail leading to a negative D^U, producing a reversal of the motion 
at (or before) R m ax- 

As mentioned above, Skripkin assumes the pressure to be isotropic. How- 
ever if we assume the evolution to be expansionfree (and fj, — /iq — constant) 
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then the isotropy of pressure follows from (80), and using (83) in (78) or (79) 
yields 

8nP = - 3k, (86) 
R 2 R K ' 

which, of course, satisfies (81). From the matching condition (22) we have 
for (86) 

SnP = e) - 3A; = 0. (87) 

R 2 R V ; 

which gives 

r 2 = kRl ie) + ci, (88) 

where c\ is an arbitrary constant. 

The mass function (18) with (44) and (83) becomes 

m = l(kR*-T 2 ) = t(R*-Rl (e) )- C j, (89) 

where we used (88). Measuring m on we obtain the total mass of the 
configuration M 

m S = M = -|. (90) 



Thus 
and 



r 2 = kR% {&) - 2M. (91) 
k 



m = -(R i -Ry ) )+M. (92) 

As mentioned before, it should be clear from physical considerations that 
the assumption of vanishing expansion (with the constant energy density 
condition) in the evolution of the fluid distribution, implies the formation of 
a vacuum cavity within the sphere. 

Applying matching conditions (23) and (24) on the boundary surface E^, 
delimiting the cavity, we obtain 

M = |(i^ (e) -i^. (i) ), (93) 

and using (93) in (91) 

T2 = kR% (i) . (94) 
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Due to the constancy of M we obtain from (93) 



2 



R*<) = ^ ^(0, (95) 
and from (44) and (95) 

A E( e)=A E(i ), (96) 

producing, because of (95), 

C/ S (e) = f|^Vc/ sW , (97) 

which implies, as expected, that the inner boundary surface E^ moves faster 
than the outer boundary surface E^. This result can also be deduced from 
the very definition of U. 

Indeed, using (44) in (27) we obtain 

u = w < 98 > 

which evaluated on E^ and E* 6 ) produces 

1~1 = C/ S (e)-R|(e) = U^R^d), (99) 

implying (97). 

Observe that from (63), (64) and (98) it follows 

° = -% ( 10 °) 

which is the solution of (77). 

It should be observed that since the pressure vanishes on EW and T^ e \ it 
should have a maximum somewhere between the two, i.e the pressure gradient 
must vanish on some spherical surface (S) within the fluid. If we denote the 
areal radius of such surface by R = Rs, then it follows from (81) 



and after integration 



{R 2 R') = (101) 



R' = f 2 (102) 
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where c 2 is a constant. Substituting (102) into (83) we have 



kR 6 + (4 -l)R A - t 2 R 3 - r\ I 0. 



(103) 



Thus the surface R$, which is the root of (103), divides the fluid into two 
regions. The inner one, with a positive pressure gradient, and the outer one 
with a negative pressure gradient. 

We are now able to prescribe the strategy to determine the Skripkin mod- 
els. First of all let us recall that without loss of generality Skripkin chooses 
T\ = R-£( e ) -R^(e) and, consequently, T\ = R^R^^ and A S ( e ) = = 1. 

Then, the integration of (83), which can only be expressed in terms of 
elliptic functions, produces 



Evaluating (85) on we obtain a differential equation for _R E(i ) whose 
integration provides its time dependence, and therefore of R(r, t) . 

If we deviate from Skripkin model and relax the condition fi = constant, 
allowing for r dependence of /i, then we need to integrate (82) instead of 
(83), which of course requires the specific r dependence of \i or m. 

Also, we could consider the anisotropic case, which allows for a t de- 
pendence of /i, in this case, of course, a specific equation of state for the 
anisotropic pressures is required (or an equivalent ansatz). 

Finally the integration in the general dissipative case would require a 
thorough knowledge of the energy production within the fluid (or a set of 
equivalent ansatze). 

7 Conclusions 

We have seen so far that expansionfree condition allows for the obtention 
of a wide range of models for the evolution of spherically symmetric self- 
gravitating systems. Ranging from nondissipative spheres, with constant en- 
ergy density and isotropic pressure (Skripkin model), to general dissipative 
anisotropic models. 

Observe that even if the Skripkin model is the simplest, from the phys- 
ical point of view, it might not be so from the mathematical point of view. 
Indeed, we could in principle choose a mass function, such that (82) could be 




(104) 
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integrated in terms of elementary functions, obviously it remains to be seen 
if such models are endowed with any physical interest. 

One of the most interesting features of the models, is the appearance of 
a vacuum cavity within the fluid distribution. It is not clear at this point 
if such models might be used to describe the formation of voids observed at 
cosmological scales (see [17] and references therein). 

The two concepts of radial velocity discussed in section 5 allows to un- 
derstand the meaning of the expansionfree evolution. As a by-product of 
such discussion, the shearfree flow (in the nondissipative case) appears to be 
equivalent to the well known homologous evolution. Particularly remarkable 
is the fact that the expansionfree ejection (collapse) implying an increase 
(decrease) in the areal radius of a layer of matter, proceeds with a decrease 
(increase) in the distance of neighbouring particles along the radial direction. 
Also, the possibility of a "splitting" of the fluid distribution (change of sign 
in U) due to dissipation and/or shear, as indicated by (74), deserves to be 
explored further. 

Finally, it is worth noticing that in the locally anisotropic case, the ex- 
pansionfree evolution, due to the second term on the left of equation (80), 
does not imply that energy density remains time independent. This situation 
becomes intelligible when it is remembered, that in the limit of hydrostatic 
equilibrium, when U = q = e = 0, we obtain from (36), 



which is just the generalization of the Tolman-Oppenheimer-Volkoff equation 
for anisotropic fluids, obtained in comoving coordinates [8]. Thus the term 
2(P r — P±) I R represents a force associated to the local anisotropy of pressure, 
and therefore the second term on the left of equation (80), is the rate of work 
done by that force, resulting in a change of /i. 
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